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ABSTRACT 


Turbulence  models  are  developed  by  supplementing  the  renormalization  group  (RNG) 
approach  of  Yakhot  &  Orszag  with  scale  expansions  for  the  Reynolds  stress  and  production 
of  dissipation  terms.  additional  expansion  parameter^^^-6?AV^)  's  the  ratio  of  the 
turbulent  to  mean  strain  time  scale.  While  low-order  expansions  appear  to  provide  an 
adequate  description  for  the  Reynolds  stress,  no  finite  truncation  of  the  expansion  for  the 
production  of  dissipation  term  in  powers  of  t]  suffices  -  terms  of  all  orders  must  be  retained. 
Based  on  these  ideas,  a  new  two-equation  model  and  Reynolds  stress  transport  model  are 
developed  for  turbulent  shear  flows.  The  models  are  tested  for  homogeneous  shear  flow 
and  flow  over  a  backward  facing  step.  Comparisons  between  the  model  predictions  and 
experimental  data  are  excellent.  F#r 
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I.  Introduction 


According  to  tlic  Kolmogorov  theory'  of  turbulence,  the  dynamics  of  velocity  fluctuations 
v{k,u,’)  at  the  scale  /  =  27r/A’  depend  only  on  the  mean  rate  of  dissipation  of  kinetic  energy 
8  and  the  length  scale  /.  Meither  the  integral  (L)  nor  dissipation  [Id  =  ‘Iir/kd)  scales  enter 
the  dynamics  of  tlie  inertial  range  of  scales  where  Id  <€.  I  L.  Based  on  this  theory,  the 
velocity  of  eddies  of  size  /  scales  as  so  that  the  characteristic  (turnover)  time  of  these 

eddies  is 

(1) 

v{l) 

The  velocity  correlation  function  can  then  be  represented  in  the  scaling  form 

<  e,(k,a;)e.(k',a;')  >~  =  C{k,io)  (2) 

\£  '  ky^J 

where  the  scaling  function  <f){x)  is  to  be  determined  Kom  other  considerations.  The  energy 
spectrum  E{k)  is  of  the  form 

E{k)  ~  k^  J  C{k,oj)dio  =  CKt^'^'k-y-^  (31 

where  CV  Bt-  so  called  Kolmogorov  constant.  These  results  also  yield  the  effective  (tur¬ 
bulent)  viscosity  uj  at  the  scale  /  =  'Iir/k: 

(4) 

This  effective  viscosity  plays  a  profound  role  in  turbulence  modelling.  For  example, 
it  is  quite  reasonable  to  assume  that  the  large-scale  properties  of  the  flow  arc  governed 
by  effective  equations  of  motion  similar  to  the  Navier-Stokes  equations  with  the  effects  of 
strong  interactions  between  the  velocity  fluctuations  taken  into  account  through  an  effective 
viscosity  (see  Yakhot  and  Orszag^) 

-Y  ~  (5) 

However,  the  strict  notion  of  eddy  viscosity  recpiires  the  existence  of  a  small  parameter 
Ij L  C  1  which  is  absent  in  turbulent  flow.  Still,  the  eddy  viscosity  concept  proves  to  be 
extremely  useful,  working  much  better  than  expected.  This  situation  is  not  unique:  in  a 
fluid  close  to  thermal  equilibrium,  the  molecular  viscosity  repre,sentation  is  very  accurate 
even  when  the  mean  free  path  A  is  not  that  small  (A/L  ~  1). 

Many  years  Ijefore  Kolmogorov’s  work,  Osborn..  Ileynclds  lealized  that  turbulent  flow  was 
r'.nrP.o  vy^ttcxTi  wiiu  li  tiaci  to  bc  treated  using  statistical  methods.  'Fbe  work  of  Reynolds 
was  similar  in  c{jncept  to  the  work  of  Boltzmann,  Gibbs  and  others  who  formulated  the 
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kinetic  theory  of  gases.  It  was  not  an  accident  that  by  the  early  20th  century  the  ideas  of 
kinetic  theory  were  adapted  to  describe  turbulence  with  turbulent  eddies  as  the  molecules 
or  building  blocks  of  this  “gas.”  By  analogy  with  kinetic  theory  the  turbulent  viscosity  is 
taken  to  be: 

~  ^rms^'  (6) 

Combining  (6)  with  (5)  and  using  €  ~  L  (so  that,  in  the  absence  of  turbulence  pro¬ 
duction,  the  turbulence  is  damped  in  one  large-eddy  turnover  time  L/urms\  we  obtain  the 
K  —  £  model: 

=  (7) 

where  K  —  |u7u7  and  is  a  constant.  The  relation  (7)  is  more  convenient  than  (6)  since  it 
expresses  the  turbulent  viscosity  through  the  directly  measurable  quantities  K  and  £. 

To  implement  (7)  we  need  equations  of  motion  for  K,S  and  the  velocity  field  v: 


dK 

dt 


9V  2 

-i-  V  •  Vv  =  -Vp  -f  i^oV  V 


-k  V  •  WK  =  -£-  V,(vip)  -f  lyoV^K 


_  +  V  •  V£  =  -Ti  -  T2  -  2:^o(Vjv0(V,V,p)  + 


where  p  is  the  pressure  and 


T,  =  2!.o(Vyu.)(V2U,)tV,t;.) 


(8) 

(9) 

(10) 


T2  = 


(11) 


Equations  (8)-(ll)  must  be  solved  subject  to  initial  and  boundary  conditions  as  well  as  the 
incompressibility  constraint  V  •  v  =  0.  In  order  to  derive  equations  for  the  mean  values  of  K 
and  £,  equations  (8)-(ll)  must  also  be  averaged  over  an  ensemble  of  the  fluctuating  part  of 
the  velocity  field  v.  In  the  absence  of  a  systematic  averaging  procedure,  turbulence  modellers 
have  used  intuitive  reasoning  combined  with  experimental  data,  tensorial  and  dimensional 
invariance,  as  well  as  scaling  arguments  based  on  the  Kolmogorov  theory^”®. 

In  the  siiiiple  ca«e  of  decaying  turbulence  with  no  mean  velocity  field,  the  equation  of 
motion  for  the  mean  kinetic  energy  is  usually  wri^t^'u  as: 


where  z/  —  i/q -f  is  the  total  viscosity  and  a/^-  is  a  proportionality  constant.  The  derivation 
of  the  equation  for  £  is  more  difficult  since  it  is  easily  shown  that  both  Tj  and  T’2  in  (11)  are 

j  y2  _ ^2  ~~ 

0(/?e.'  )  where  Rc,  =  K  /uo£  =  0{  Urms^li^o)  is  the  turbulence  Reynolds  number.  However, 
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the  following  bold  assumption  (cf.  Tennekes  and  Lumley^)  is  fruitful;  it  is  assumed  that  the 

singular  contributions  to  Tj  and  cancel  each  other  in  the  limit  as  Rtt  — >  oo  and  the 

remaining  0(1)  part  can  be  written  as 

T^+n  =  Cs2=  (13) 

where  0^2  =  0(1).  Thus,  in  decaying  turbulence. 

dS  ^  s'"  d  (  d£\ 

dt  ^  dxi  ^  ^ 

where  as  =  0(1).  The  assumption  (13)  was  not  rigorously  justified  but  it  has  led  to  simple 
equations  useful  for  practical  calculations.  In  the  case  of  homogeneous  isotropic  turbulence 
the  diffusion  terms  in  (12)  and  (14)  disappear  yielding 

K  oc  (1-5) 

where  7  =  i/(C^'2  —  1).  The  fact  that  the  coefficient  Cs2  determines  the  power  law  of 
turbulence  decay  demonstrates  the  significant  role  of  dimensionless  constants  in  the  theory 
of  turbulence.  Indeed,  even  a  small  error  in  Ce2  is  amplified  in  the  limit  as  f  >  00.  The  first 
term  on  the  right  side  of  (14)  has  a  simple  interpretation  in  terms  of  a  relaxation  time: 

c  S 


where 


e 


is  the  only  turbulence  time-scale  that  one  can  construct  from  the  parameters  of  the  problem. 

In  flows  with  non-zero  values  of  the  mean  strain  5, 7  =  |(|^  +  |^),  the  velocity  field  in 
(8)-(ll)  can  be  decomposed  into  mean  (U)  and  fluctuating  (u)  components.  In  this  case  the 
equations  of  motion  are  typically  taken  to  be  of  the  form: 


— -bU-VA'  =  -ro5'.,-<f+—  aKi' 


d£  ^  e  ^  ^  e\  d  (  d£\ 

-  +  U  .  VE  =  -Cn^r„s.,  -  C„=  +  —  (“s-g;-) 

wL:--  =  IZfuT  is  the  Reynolds  stress  tensor  and  the  coefficients  Ce\  and  (7^2  are  constants 

that  are  usuallv  d<’tciinincu  fiou'  benchiiianv  ^.xp  crimen  Is.  the  t  -  production  term  in  (16) 
comes  from  the  contributions  to  the  S  -  equation  (10)  of  the  type; 

^  du,  duj  ^ 

Q  r'.  aij 

axi  0x1 
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which  are  modeled  as 


in  the  relaxation  time  approximation. 

There  are  also  d3’aamical  constraints  that  a  consistent  turbulence  model  must  obey.  First 
of  all,  it  is  obvious  that 

Z>0,  T>0.  (17) 

It  can  be  shown  that  the  conditions  of  realizability  (17)  are  satisfied  for  homogeneous  tur¬ 


bulent  flows  if: 


when 


K,e  0. 


Secondly,  the  turbulence  model  must  be  invariant  under  the  Galilean  transformation: 

X*  =  X  —  Uo<  (18) 

where  Uq  is  any  constant  velocity.  In  this  sense,  the  turbulence  model  must  behave  like 
the  original  Navier-Stokes  equations  which  are  Galilean  invariant.  Any  model  violating 
invariance  under  (18)  is  physically  incorrect. 

Using  the  simplest  closure  for  the  deviatoric  part  of  r,j: 

Tij  =  ~2uTSij  =  -2C^^Sij  (19) 

(where  ^  0.085  is  a  constant)  (16)  is  transformed  into  a  familiar  form.  It  can  be  shown 
that  (16)  with  (19)  satisfies  the  realizability  conditions  /f  >  0  and  T  >  0  provided  that 
Ce2  ^  1  and  >  0.  These  and  other  interesting  properties  of  equations  (16)  and  (19)  have 
been  recently  reviewed  by  Speziale®. 

Let  us  introduce  the  dimensionless  variables 

SK  K 

where  S  —  (2S'ij5’p)‘^^  and  Kq  is  the  initial  turbulent  kinetic  energy.  In  homogeneous  shear 
flow  where  V^K  =  =  0,  equations  (16)  and  (19)  have  the  simple  form: 

-  h  -  Tf  ) 


^  =  -C,{Cs,~l)v^  +  {Ce2-l) 
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where  t'  —  St.  Equation  (20)  has  a  single  fixed  point  given  by 


which  is  obtained  by  setting  drj/dt  =  0.  This  fixed  point  yields  asymptotic  solutions  for  K 
and  S  of  farm 


K  (X  e' 


oc 


where  the  growth  rate  A  is  given  by® 


(ai-l)(Ce2-l) 


Physical  and  numerical  experiments  on  homogeneous  shear  flow  -  for  not  too  large  values  of 
the  dimensionless  shear  rate  r]  -  indicate  that  indeed,  rj  rjo  and  K,£  oc  e"'*’  for  t*  '>  1. 
This  means  that  the  solutions  of  any  turbulence  model  must  be  attracted  to  this  fixed  point 
when  Tj  is  not  too  large.  The  fact  that  the  simple  K  —  £  model  is  capable  of  describing  such 
a  non-trivial  behavior  is  remarkable. 

Another  important  consequence  of  equations  (16)  and  (19)  is  the  Reynolds  number  inde¬ 
pendence  of  the  von  Karman  constant  in  turbulent  channel  flow.  It  is  easy  to  show  that  the 
normalized  dimensionless  velocity  profile  U'^  is  given  in  the  region  of  constant  energy  K'^ 
by; 

=  Un  B  (23) 

K 

where  k  is  the  von  Karman  constant.  The  dimensionless  variables  ?/+  and  U'^  will  be  discussed 
in  more  detail  in  Section  IV.  It  follows  from  the  model  (16)  and  (19)  that  the  von  Karman 
constant  is  independent  of  Re.  This  is  the  result  of  the  cancellation  of  the  singular  terms 
assumed  in  the  derivation  of  (13).  If,  on  the  other  hand  +  T2  =  O(Rey^)  and  the 
assumption  (13)  is  incorrect,  then  the  von  Karman  constant  k  and  the  constant  B  in  (23) 
will  depend  on  Ret.  This  would  have  a  significant  impact  on  engineering  calculations  since 
(23)  enters  in  the  expression  for  the  friction  coelRcient.  Surprisingly,  even  today,  there  is 
no  compelling  experimental  evidence  for  the  constancy  of  B  and  k,  although  the  prevailing 
opinion  is  that  B  and  k  are  constants. 

Now  let  us  discuss  the  closure  assumption  (19),  which  can  be  rewritten  as; 

-  -2C,J\7/,j  (24) 

where  =  Si^K !£  is  a  dimensionle.ss  tensor  which  is  equal  to  zero  in  isotropic  turbulence. 
When  rjij  is  small,  the  deviation  from  the  isotropic  solution  is  small.  Unfortunately,  in  many 


.5 


practical  applications  of  the  K  —  £  model,  like  channel  flow,  the  values  of  parameter  t/,j  vary 
from  //  ~  20  in  the  buffer  region  to  r/  ~  3  in  the  logarithmic  layer  (where  rj  =  llr/ijlD- 

The  dimensionless  tensor  can  be  used  for  a  perturbation  expansion  of  fij  in  powers 
of  7;  when  the  departures  from  isotropy  are  small; 

r,j  =  ~K[2C^r^,j  +  (9(7/5)  +  ...].  (25) 

In  this  case,  the  theory  of  strained  turbulence  can  be  formulated  in  terms  of  a  double  expan¬ 
sion  in  powers  of  the  two  dimensionless  small  parameters  Re~^  and  T/p  .  This  expansion  gives 
rise  to  the  anisotropic  eddy  viscosity  explored  by  Speziale®  and  Rubenstein  and  Barton'°. 
The  existence  of  the  second  dimensionless  expansion  parameter  7/  is  not  reflected  in  the 
functional  form  of  the  standard  K  —  £  model.  This  may  account  for  the  relatively  poor 
performance  of  the  standard  model  in  flows  with  large  values  of  r/p. 

It  is  apparent  that  a  systematic  derivation  of  turbulence  models  will  enable  us  to  address 
some  of  following  questions: 

(i)  Is  it  true  that  Ti  +  T2  =  (9(1)  as  in  (13)? 

(ii)  Does  the  model  (16)  include  all  of  the  relevan!  effects  or  is  something  missing? 

(iii)  Is  it  possible  to  derive  a  more  complete  Galilean  invariant  model  than  (16)  satisfying 
the  constraints  discussed  above,  namely: 

(a)  A'  >  0;£'  >  0  (realizability) 

(b)  dq/dt  —  0  when  7/  -4  7/0  (fixed  point) 

(C)  K  =  0(1)? 

Condition  3(a)  and  the  reciuirement  of  Galilean  invariance  are  the  basic  constraints,  while 
constraints  3(b)-(c)  are  based  on  experimental  data  and  must  be  accepted  with  some  measure 
of  skepticism  until  theoretically  justified. 

II.  Renormalization  Group  Methods  and  Turbulence  Models 

The  equations  of  motion  (8)-(ll)  describe  general  turbulent  flow.  We  assume  that  the 
turbulence  is  driven  by  mean  strains  and  decompose  the  velocity  field  into  mean  and  fluc¬ 
tuating  parts  (v  =  U  +  u)  as  in  Section  I.  We  also  make  the  basic  assumption  that  the 
turbulence  statistics  are  homogeneous  when  the  mean  strain  .Sp  is  not  too  large.  Next,  we 
assume  that  the  fluctuating  velocity  field  u  is  governed  by  the  Navier-Stokes  equations  driven 
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by  a  random  force,  chosen  in  such  a  way  that  the  global  properties  of  the  resulting  held  are 
the  same  as  those  in  the  flow  driven  by  the  m^an  strain  S'  . 

To  derive  the  equations  for  u,  we  consider  an  infinite  domain  in  which  a  Newtonian  fluid 
is  stirred  by  a  Gaussian  random  force.  The  force  is  defined  by  its  correlation  function  in 
wavevector  and  frequency  space  (see  Yakhot  and  Orszag'^), 


2Do{27ry+^k-yp,^{k}6{k  +  k')6{u;  +  a-'),  <  k  <  oc 

0 ,  A  5  <  A'  <  A  L 


(26) 


where  k  is  the  wavevector,  u  is  the  frequency,  k  =  jkj  and  d  is  the  number  of  space  di¬ 
mensions.  The  parameter  y  is  chosen  to  give  the  Kolmogorov  form  of  the  small  scale  en¬ 
ergy  spectrum;  and  in  three  dimensions  y  =  d  =  3.  Statistical  homogeneity  in  space  and 
time  is  guaranteed  through  the  factors  6(k  +  k')(S(u;  -f  u)')  in  (26).  The  projection  operator 
Pij{k)  —  —  kikj/P  makes  the  force  statistically  isotropic  and  divergence-free.  In  the  limit 

of  infinite  Reynolds  numbers,  i\l  ->  0  (but  A^,  >>  A,,  ^  0).  The  forcing  function  (26)  re¬ 
flects  the  fact  that  turbulence  is  usually  driven  by  hydrodynamic  instabilities  with  the  most 
unstable  mode  at  k  =  Al.  The  traditional  view  is  that  this  primary  instability  generates 
velocity  fluctuations  with  k  >  \l  such  that  there  is  a  direct  energy  cascade  to  small  scales. 
The  energy  cascade  may  itself  generate  a  universal  random  force  of  the  form  (26),  but  this 
has  not  yet  been  demonstrated  rigorously. 

Our  model  is  the  forced  Navier-Stokes  equations: 

dvi  _  _2  /• 

~  +  VjVjV,  =  -V.p  -t-  +  /, 

V,c,  =  0  (27) 


where  the  Gaussian  force  f  is  given  by  (26)  and  the  density  p  has  been  absorbed  in  the 
pressure  p.  In  contrast  to  th.  earlier  development  of  RNG  models  in  Ref.  (2),  the  only  new 
feature  of  the  model  (26)-(27)  is  the  infrared  cutoff  of  the  random  force:  <  fifj  >=  0  when 
0  <  A’  <  A/,.  This  property,  which  is  usually  unimportant,  is  needed  if  we  wish  to  derive  an 
equation  foi  the  mean  rate  of  energy  dissipation  S . 

The  dynamical  equations  for  the  total  kinetic  energy  per  unit  mass  K  =  and  the 

homogeneous  part  of  the  instantaneous  rate  of  energy  dissipation  per  unit  mass  £  = 
are  obtained  from  (27), 


(28) 

(29) 
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In  the  equation  for  £,  driven  by  large-scale  features  such  as  boundary  and  initial  con¬ 
ditions,  1\  balances  T2  to  leading  order.  .\s  we  shall  see.  for  st<'ady-state  flows  (27)-(29)  is 
sustained  by  the  torce  f;  7’.  is  balanced  by  both  /’i  and  ‘he  random-force  contribution  to  the 
if-production  given  by  7^-  =  'Iv'ol'V 

We  seek  equations  tor  the  mean  values  U  =<  v  >.  K  =<  K  >  and  £  =<  £  >,  averaged 
over  an  ensemble  of  the  fluctuating  velocity  field.  lo  find  these  equations,  we  shall  use  the 
dynamic  renorm.alization  group  and  the  £-e.\pansion.  Since  the  renormalized  eciuations  for 
A  and  £  may  not  be  trivially  related  to  the  renormalized  .\a\’ier-Stokes  equations,  the  HNG 
procedure  must  be  applied  to  all  of  the  eriualions  (27)-(29). 

The  renornic'ilizatio.']  group  appli('d  to  (28)-(29)  is  de.scribed  in  detail  in  ^’akhot  and 
Smith”.  IIer<'  we  present  the  main  results.  I  he  equations  of  motion  averaged  over  the 
fluctuating  velocity  field  v  are; 


—  +  U  .  VA  =  -  r„,s,,  +  — 


as  I'-: 


()£  £  __  ,,  £^  _  77  /  d£ 

Ot  A  A  rdx,  \  d.;q 


^  duidui  ^ 

7^  =  2/zo8q  — 7] — ,  (31) 

ax,  dxj 

dk'i  =■  1.‘12,  Ckq  =  1.68,  and  ok  =  cvf  =  1.39.  The  relations  (30)-(31)  are  obtained  in  the 
lowest  order  of  the  effective  Reynolds  number: 

^  oc£~0(l)  (32) 


where  .A,  —  ‘lx j L  and  v  is  the  effective  viscosity  (7).  The  parameter  e  =  7  — y  is  characteristic 
of  the  random  noise  driving  this  artificial  flow.  It  was  also  shown”  that  in  this  order  of 
the  expansion,  T\  -f-  T2  =  0(1)  which  confirms  the  hypothesis  leading  to  (13).  In  sheared 
turbulence  where  ^  0,  the  Reynolds  decomposition  generates  contributions  of  the  order 

e  77iqc»u/ 
dx  ' 

\J  mL  ^  \J  Ji-  J 


d’he  isotropic  and  homogeneous  random  velocity  field  v  was  assumed  to  be  governed 
by  the  equations  (2)-(4).  The  renormalization  group  and  the  t-expansion  were  applied 
1,0  evaluate  (33).  'fhis  .systemmatic  procedure  generated  an  expansion  in  powers  of  the 
dimensionless  parameter  r/q  =  S,jl\  l£.  The  contribution  to  TZ  in  (31)  is  small  for  weakly 
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siraiut'd  turbulence  and  large  in  the  rapid  distortion  limit  when  //  — ■  oc.  dil  attempts  to 
cdose  (31;  using  the  methods  based  on  t he  c-expansion  were  fruitless.  In  tin'  next  s('clion  we 
shall  show  how  to  evaluate  (31)  by  applying  tin-  geuieral  criteria  outlined  in  Section  1. 


III.  Generalized  A  —  £  model 


.\eg!ecting  'K  in  (31)  and  using  the  low  order  closure  (19)  for  we  obtain  equations  of 
motion,  which  are  those  of  the  standard  A  - 

Dl<  - 


Dt 


D£ 

IJt 


£_ 

1< 


‘j 


model: 

d  ( 

OK 

■  ay 

dx^ 

T  ± 

/ 

"  K  ^  d.r, 

1“-' 

where  D/  Dt  =  O/Ot  +  U  •  V. 

This  model  satisfies  all  of  the  criteria  formnlated  in  Section  I.  In  homogeneous  shear  flow 
wheie  (j''Jdxj  =  SS^ySyi.  eciuations  (34)  have  the  fixed  point  r;o  =  [c ^“"-1)]  ^ 
for  St  >  1  the  kinetic  energy  and  dissipation  grow  exponentially: 


Ao  ' 

with  a  growth  rate  ^  =  0.112,  wh’ch  is  very  close  to  tiiat  obtained  from  numerical  and 
ohysical  experiments.  These  results  wdll  be  discussed  in  more  detail  in  Section  V. 

The  model  (31)  is  based  on  the  assumption  that  //  is  small;  however,  the  numerical  value 
of  //  1.38  does  not  satisfy  tliis  constraint.  Furthermore,  in  the  logarithmic  region  of  a 

turbulent  channel  flow',  the  value  of  the  dimensionless  shear  ?/  =  SK /£  'Uc  3  is  also  not  small. 
Ihis  mc’ans  that  the  contribution  uo  Eq.  (30)  due  to  TZ  cannot  be  neglected  there. 

Iterating  the  expression  for  TZ  using  the  Navier-Stokes  eciuations  w'ill  generate  a  power 
series:  _ 

K  =  >'T.eE.'.(^')'  US) 

I!=0  \  / 

vvtrere  S  —  (2.S,j .S',^)^^,  It  is  clear  that  w'hen  .S'  — ^  0,  7^  goes  to  zero  faster  than  and 

may  be  neglected  in  (30).  To  use  (.35)  in  llie  case  of  large  shear,  let  us  consider  as  a  model 
only  the  subset  of  terms  corresponding  to  the  geometric  series 


7^"  = 


T(' 

T|=0 


3n 
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c 

c 


(36) 


We  do  not  know  the  values  of  the  coefficients  aiui  cannot  evaluate  the  magnitude  of 
R  —  7v,\  flowever,  this  form  of  is  plausahle  since  (36)  implies  that 


O 


for  all  values  o‘  //.  riius.  'R^  vauislu's  when  <f  0,  a  constraint  that  also  follows  from 
the  clelinitioii  of  R  given  in  (31),  so  that  {3u)  maintains  the  satisfaction  of  the  realizability 
conditions  (17).  Assuming  tliat  tiie  fixed  point  value  //o  —  4.38,  obtained  in  the  limit  of 
small  .Wj.  is  inva''iant  to  dropping  all  term.s  bu^  th  in  (36)  we  postulate  that: 

lyrS-il  - 

R  =  - ^—7^.  (37) 

1  +  Jr 

Fajuatiou  (37)  can  be  rearranged  into  the  alternative  form 

c  r  ("i  -  i) 

R  =  1  7  T 

1  +  Jtj-'  k 

from  which  it  i.s  clear  that  the  crucial  constraint  that  'R  — ►  0  as  £1  — >•  0  is  satisfied.  The 
limits  of  w(\,k  and  .strong  strains  can  also  be  easily  discerned  from  (38).  It  is  clear  that  in 
the  limit  of  weak  strains  where  i]  — >■  0: 


0- 

h 

On  the  other  hand,  in  the  limit  of  strong  strains  where  i]  — >  oc: 

3  r  R 

which  b('comcs  singu'ar.  From  the  definition  of  R  given  in  (31),  it  follows  that 


where' 


^  dJCt 


It  is  a  simple  matter  to  show  Tiat  ||.4,j1|2  <  |  and  —  2  limit  as  t]  — >  cc, 

it  hdiows  diat  R  =  (Mjl)  consistent  with  (46)  Furthermore,  in  the  limit  as  o  — >  0,  7^  — >  0 
as  0{£‘']  and  in  the  limit  as  t/  — >  0,  "R.  0  a.s  Oir/^)  or  0[if  )  (which  follows  from  local 

i.^otropy;  c.f.  .Monin  and  Yaglom''^).  4  hese  results  are  generally  consistent  with  (39). 

It  i.s  clear  that  the  model  (38),  while  not  rigoromly  derived  by  RNG  methods,  satisfies 
a  variety  of  important  consistency  conditions  tliat  [ilace  .severe  constraints  on  the  possible 
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functional  form  of  7^.  In  particular,  no  finite  number  of  terms  in  the  series  (35)  can  give  the 
desired  asymptotics.  While  the  analyticity  of  TZ{7])  in  the  limit  as  r/  oo  is  not  necessary 
we  choose  to  assume  it  for  sim.plicity.  If  this  is  the  case,  then  the  postulate  (38)  for  TZ  is 
quite  reasonable. 


IV.  The  von  Karman  Constant 

The  proposed  model  for  TZ  given  in  (38)  contains  one  undetermined  constant  jS  which  will 
now  be  related  to  the  von  Karman  constant.  In  order  to  do  this,  we  rewrite  the  dissipation 
rate  transport  equation  in  the  equivalent  form 

f  +  U  .  VJ  =  -CS.L..A, 


vhere  the  coeflicient  C*2  given  by 


(l  -  f ) 

=  1.68+ - .  . 


1+,V 

In  the  log  layer  of  a  turbulent  boundary  layer  the  turbulence  production  — r,j5jj  =  S  and, 
hence,  rj  =  1/^^  which  renders  CI2  constant.  Furthermore,  convection  effects  (i.e.,  DS/Dt) 
can  be  neglected  and  the  mean  velocity  assumes  the  logarithmic  form  (23)  where 

ly  Ur 

and  Ur  is  the  shear  velocity.  With  these  assumptions,  it  is  a  straightforward  calculation  to 
show  that  (42)  yields  the  constraint'^ 

K  =  - - -  (44) 

Of  I 

v.'hich  establishes  a  direct  relationship  between  the  von  Karman  constant  a  and  the  constant 
/?  in  the  model  (38)  for  TZ.  It  follows  that  when 


13  =  0.012 


the  von  Karman  eonstant  k  =  0.4.  Fortunately,  k  is  not  very  sensitive  to  variations  in  /3  (in 
applications,  can  be  varied  between  0.01  and  0.015  with  only  small  changes  in  k).  With 
this  selection  of  /i  (and  the  fixed  point  7/0  =  4.38)  we  now  have  a  closed  model  for  TZ.  Since 
the  model  yields  a  log  layer  with  a  von  Karman  constant  k  of  0.4,  it  will  give  good  results 
for  turbulent  channel  flow.  The  model  will  be  subjected  to  a  more  severe  wall-bounded  flow 
test  in  Section  VI  where  turbulent  flow  over  a  backward  facing  step  will  be  considered. 
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V.  Homogeneous  Shear  Flow  and  the  Relaxation  Time  Approximation 


In  this  section  we  report  tests  of  the  new  model  in  homogeneous  shear  flow  where  an 
initially  isotropic  turbulence  (with  turbulent  kinetic  energy  Ko  and  dissipation  rate  Fo)  is, 
at  time  t  =  0,  subjected  to  a  constant  shear  rate  5  with  the  corresponding  mean  velocity 
gradient  tensor 

— -  —  bdiidj2. 

dxj 

The  model  given  by  equations  (19),  (30),  and  (38)  yields  a  simple  set  of  coupled  nonlinear 
ordinary  differential  equations  since  K  =  K{t)  and  £  =  £{i).  The  solution  obtained  for  the 
time  evolution  of  the  turbulent  kinetic  energy  is  compared  in  Figure  1  with  the  large-eddy 
simulation  of  Bardina  et  al.^‘‘  Despite  the  fact  that  the  model  yields  an  excellent  long-time 
growth  rate  for  K  (i.e.,  the  model  predicts  that  A  —  0.142  in  comparison  to  physical  and 
numerical  experiments'^''®  which  indicate  that  A  is  in  the  range  of  0.14  -  0.16),  the  early  time 
values  of  K  are  overpredicted.  This  results  from  the  use  of  an  eddy  viscosity  representation 
for  fij  which  responds  too  quickly  to  the  application  of  the  shear  at  t  =  0. 

The  early  time  behavior  can  be  described  more  accurately  when  a  relaxation  time  ap¬ 
proximation  is  introduced.  In  this  approximation,  we  allow  the  Reynolds  stress  r,j  to  relax 
to  the  eddy  viscosity  model  (19)  as  follows: 

d  J 

—{Tij  -f  2i/r5q)  =  -Ct„=(D;  +  2uTSij)  (45) 

for  homogeneous  flows  where  Cth  is  a  relaxation  coefficient.  It  is  a  simple  matter  to  show 
that  as  a  direct  consequence  of  (45) 


^  +  C,2  -  2)KS„  +  0(5J).  (46) 

When  an  initially  isotropic  turbulence  is  strained.  Crow"’  showed  that 

CACt,  +  O,  -  2)  =  4 

which  yields  ~  3.5.  Eq.  (46)  gives  a  simple  Reynolds  stress  transport  model. 

In  addition  to  (45)  for  the  relaxation  time  approximation,  should  be  replaced  by 
the  production  P  =  in  Eq.  (37)  giving  the  relaxation  model 


(47) 


where  for  homogeneous  shear  flow  P  =  — ri2.S’.  When  the  eddy  viscosity  model  (in  which 
P  =  2urSij£tj)  is  introduced  in  (47),  Eq.  (37)  is  recovered. 
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In  Figure  2,  the  time  evolution  of  the  turbulent  kinetic  energy  in  homogeneous  shear 
flow  obtained  from  the  relaxation  model  (46)  -  (47)  is  shown.  These  results  are  in  excellent 
agreement  with  the  large-eddy  simulation  of  Bardina  et  al.^'^  This  is  quite  encouraging  con¬ 
sidering  the  fact  that  this  relaxation  model  was  calibrated  independent  of  homogenous  shear 
flow  (i.e.,  it  was  calibrated  by  setting  the  value  of  the  von  Karman  constant  to  0.4  and  by 
using  the  Crow'  result  for  plane  strain  turbulence). 


VI.  Turbulent  Flow  over  a  Backward  Facing  Step 


Turbulent  flow  over  a  backward  facing  step  is  one  of  the  standard  test  cases  used  to 
evaluate  the  performance  of  turbulence  models  in  separated  flows.  The  Kim,  Kline  and 
Johnston'*  experimental  test  case  will  be  considered  here  for  the  backstop  problem.  For  this 
flow  configuration,  the  expansion  ratio  (step  height:  outlet  channel  height)  F  is  1:3  and  the 
Reynolds  number  Re  =  132,000  based  on  the  inlet  centerline  mean  velocity  and  the  outlet 
channel  height  (see  Figure  3).  We  will  use  the  data  for  this  test  case  as  updated  by  Eaton 
and  Johnston'®  for  the  1980-81  AFOSR-HTTM  Stanford  Conference  on  Turbulence.  Since 
the  turbulence  is  statistically  steady,  we  can  use  the  simpler  form  of  the  model  given  by  (37) 
(i.e.,  there  is  no  need  to  incorporate  relaxational  effects).  The  mean  turbulence  equations 
are  solved  numerically  by  a  finite  volume  method  on  a  200  x  100  variable  mesh  (see  Thangam 
and  Speziale^®  for  more  details  on  the  numerical  method).  Inlet  profiles  for  U,K,  and  £  are 
specified  five  step  heights  upstream  of  the  step  corner  (a  developing  turbulent  channel  flow 
calculation  is  carried  out  that  matches  the  experimental  inlet  mean  velocity  profiles  of  Eaton 
and  Johnston'®);  conservative  extrapolated  outflow  conditions  are  applied  thirty  step  heights 
downstream  of  the  step  corner.  Along  the  upper  and  lower  walls,  the  law  of  the  wall  (23)  with 
K  =  0.4  and  F  =  5  is  used  along  with  standard  log  layer  boundary  conditions  for  the  kinetic 
energy  and  dissipation  obtained  from  a  production-equals-dissipation  equilibrium  hypothesis 
and  the  constraint  of  vanishing  normal  derivative  of  K  at  the  wall.  It  was  demonstrated 
recently  by  Thangam  and  Speziale^®  that,  provided  the  flow  field  is  well  resolved,  the  use  of 
these  law  of  the  wall  boundary  conditions  leads  to  only  minimal  errors. 

Two  sets  of  results  will  be  presented  based  on  the  algebraic  model  (37)  for  TZ:  one  set 
obtained  by  using  the  isotropic  eddy  viscosity  model  (19)  and  another  obtained  from  the 
anisotropic  eddy  viscosity  model  of  Speziale®  given  by 


-2C^-^S,j  -  ACdCI~{S,j  --  hk 
c  E 

-\-SikSkj  —  -SkiSkiSij) 


(48) 
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where  Sij  the  Oldroyd  derivative  of  Sij  with  convective  effects  neglected;  Co  is  a  dimen¬ 
sionless  constant  which  assumes  a  value  of  1.68.  For  thin  shear  layers,  this  model  is  very 
similar  to  the  anisotropic  eddy  viscosity  model  derived  by  Rubenstein  and  Barton^”  bcLsed 
on  RNG  methods.  This  type  of  anisotropic  eddy  viscosity  model  is  obtained  when  terms  of 
order  are  maintained  in  the  expansion  (25)  for  T,j. 

In  Figures  4(a)-(b),  the  computed  mean  velocity  streamlines  and  mean  velocity  profiles 
obtained  using  the  isotropic  eddy  viscosity  model  (19)  are  compared  wdth  the  experimental 
data.^®  Reattachment  is  predicted  at  xafH  ^  6.6,  a  result  that  is  within  6%  of  the  ex¬ 
perimentally  measured  reattachment  point  xnf  H  ^  7.0.  The  computed  turbulence  intensity 
and  turbulence  shear  stress  profiles  are  compared  with  experimental  data  in  Figures  5(a)-(b). 
The  agreement  is  comparably  good. 

By  combining  the  model  (37)  for  7Z  wdth  the  anisotropic  eddy  viscosity  model  (48),  even 
better  results  are  obtained.  In  Figures  6(a)'(b),  the  predicted  mean  velocity  streamlines 
and  mean  velocity  profiles  are  shown  to  compare  extremely  favorably  with  the  experimental 
data.^^  The  model  predicts  reattachment  at  xr/H  %  7.0  which  is  essentially  the  same  as 
the  experimental  result.  Furthermore,  the  model  predicts  a  noticeable  secondary  separation 
bubble  below  the  corner  of  the  step  consistent  with  experimental  observations  for  this  back- 
step  flow.  The  agreement  between  the  model  predictions  and  the  experimental  data  for  the 
turbulence  intensity  and  turbulence  shear  stress  profiles  is  also  excellent  as  shown  in  Figures 
7(a)-(b).  The  quality  of  these  predictions  is  quite  remarkable  for  a  two-equation  model. 

VII.  Concluding  Remarks 

The  renormalization  group  formalism  of  Yakhot  &  Orszag  for  the  development  of  turbu¬ 
lence  models  has  been  supplemented  with  scale  expansions  for  the  Reynolds  stress  and  pro¬ 
duction  of  dissipation  terms.  Here,  the  extra  expansion  parameter  is  taken  to  he  rj  =  SK/6 
which  is  ratio  of  the  turbulent  to  mean  strain  time  scale.  For  the  Reynolds  stress,  this 
approach  -  which  leads  to  the  development  of  anisotropic  eddy  viscosity  models  as  well 
as  Reynolds  stress  transport  models  -  is  analogous  to  that  introduced  by  Rubenstein  and 
Barton. However,  the  present  method  is  completely  new  for  the  modelling  of  the  production 
of  dissipation  term  TZ  which  is  neglected  in  most  of  the  commonly  used  turbulence  models. 
The  interesting  result  for  TZ  is  that  no  finite  order  truncation  of  the  expansion  satisfies  the 
necessary  physical  constraints;  terms  of  all  orders  in  the  expansion  must  be  retained.  This 
complication  eliminates  the  possibility  of  determining  the  model  explicitly  in  closed  form. 
However,  a  highly  plausible  form,  with  only  one  undetermined  constant,  is  postulated  here 
which  satisfies  all  of  the  necessary  physical  constraints  (i.e.,  realizability  and  consistency  with 
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the  weak  and  strong  strain  limits).  The  constant  is  calibrated  by  setting  the  von  Karman 
constant  k  to  0.4. 

The  new  models  have  been  tested  for  homogeneous  shear  flow  and  for  flow  over  a  back¬ 
ward  facing  step.  E.xcellent  results  are  obtained  in  both  cases.  For  the  case  of  homogeneous 
shear  flow,  the  best  results  are  obtained  from  the  Reynolds  stress  transport  model  (i.e.,  the 
relaxational  model  discussed  in  Section  V).  On  the  o'^^er  hand,  excellent  results  are  ob¬ 
tained  with  eddy  viscosity  models  for  the  backstop  problem.  In  all  of  these  calculations,  no 
ad  hoc  adjustments  of  the  constants  are  made.  The  applications  considered  in  the  paper  are 
restricted  to  simple  shear  flows  since  the  current  version  of  the  dissipation  rate  transport 
equation  has  only  been  modelled  to  account  for  the  effects  of  irrotational  strains.  Incorpora¬ 
tion  of  the  effects  of  rotational  strains,  which  can  be  important  in  turbulent  flow's  involving 
curvature  or  a  system  rotation,  is  a  difficult  task  that  has  not  yet  been  achieved.  The  reduc¬ 
tion  in  the  energy  cascade  that  occurs  in  rotating  isotropic  turbulence  and  the  stabilizing 
or  destabilizing  effects  of  rotations  on  homogeneous  shear  flow  are  but  two  examples.^* 
This  more  difficult  problem  of  accounting  for  rotational  strains  using  a  comparable  scale 
expansion  technique  will  be  the  subject  of  a  future  study. 
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Turbulent  flow  over  a  backward  facing  step 
Figure  3 
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BACKWARD-FACING  STEP: 


(b)  Dimensionless  mean  velocity  profile 

( - Computations  with  isotropic  eddy  viscosity; 

0  Experiments  of  Kim  et  af,  1980;  Eaton  &  Johnston,  1981) 


Computed  mean  flowfleld  for  the  new  RNG  K-e  model 
[E  =  1:3;  Re  =  132,000;  200x100  mesh] 

Figure  4 
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BACKWARD-FACING  STEP: 


(a)  Turbulence  Intensity 
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(b)  Turbulence  shear  stress 

Computed  turbulence  stressesforthe  new  RNG  K-e  model  [E = 1 :3;  Re  =  1 32,000; 

200x100  mesh;  -  computations  with  isotropic  eddy  viscosity; 

0  experiments  of  Kim  etal,  1980;  Eaton  &  Johnston,  1981] 

Figure  5 
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BACKWARD-FACING  STEP: 


X/H 

(a)  Streamlines 
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(b)  Dimensionless  mean  velocity  profile 

( - Computations  with  anisotropic  eddy  viscosity; 

o  Experiments  J  Kim  eta!,  1980;  Eaton  &  Johnston,  1981) 


Computed  mea.*  flowfield  for  the  new  RNG  K-t  model 
[E  =  1:3;  Re  =  132,000;  200x100  mesh] 

Figure  6 
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BACKWARD-FACING  STEP: 


Y 

H 


0.0  0.1  0.2  0.1  0.2  0.1  0.2  0.1  0.2  0.1  0.2  0.1  0.2 


--  1/2 

ruu\ 


(a)  Turbulence  intensity 


Y 

H 


0.0  0.01  0.02  O.ni  0.02 


0.01  0.02  0.01  0.02 


10.33  15.66 


r 

) 

) 

o 

D 

-  ^ 
o\ 

1 

/  1  - 1  _r  1 ,  1 _ 

0.01  0.02  0.01  0.02 


(b)  Turbulence  shear  stress 

ComputedturbulentestressesforthenewRNG/C-Emodel[E=1:3;Re  =  132,000; 

200x100  mesh;  -  computations  with  anisotropic  eddy  viscosity; 

0  experiments  of  Kim  et  al,  1980;  Eaton  &  Johnston,  1981] 

Figure  7 
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